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Abstract 



^ . The role of superselection rules for the derivation of classical probability 

^ ' within quantum mechanics is investigated and examples of superselection 

rules induced by the environment are discussed. 



1 Introduction 

One of the puzzles of quantum mechanics is the question, how classical objects 
can arise in quantum theory. Quantum mechanics is a statistical theory, but 
its statistics differs on a fundamental level from the statistics of classical objects. 
The EPR problem and the violation of Bell's inequalities are consequences of this 
fact [1-3] (English translation of [2] in [4]). These questions are usually discussed 
as problems of the interpretation of the Schrodinger wave function. Here I would 
like to emphasize a more abstract approach which relates these problems to the 
geometrical structure of the quantum mechanical state space, i.e. the total set 
of pure and mixed states. It is exactly this structure which forbids to assign an 
objective probability measure to the states. 
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It is known since a long time that the statistical results of quantum mechanics 
become consistent with a classical statistics of "facts", if the superposition prin- 
ciple is reduced to "superselection sectors", i.e. coherent orthogonal subspaces of 
the full Hilbert space. The mathematical structure of quantum mechanics and of 
quantum field theory provides us with only a few "superselection rules" , the most 
important being the charge superselection rule related to gauge invariance, see 
e.g. [5,6] and the references given therein. But there are definitively not enough 
of these superselection rules to understand classical properties in quantum theory. 
A possible solution of this problem is the emergence of effective superselection 
rules due to decoherence caused by the interaction with the environment. These 
investigations - often related to a discussion of the process of measurement - 
started in the eighties; some early references are [7-9]. 

In this talk the transition from quantum mechanics to classical physics is in- 
vestigated in the following restricted sense: How can a classical statistical theory 
emerge from quantum probability? The emphasis is laid on the principles of in- 
duced superselection. Many aspects, like Markov approximation and localization 
models arc totally omitted. For a comprehensive investigation of decoherence see 
the book [10]. The talk is organized as follows. The structure of the state space 
is discussed in Sect. 2. Superselection rules are defined in Sect. 3. In Sect. 4 the 
dynamics of subsystems and the emergence of effective superselection rules are 
recapitulated. Some exactly solvable models and mathematical formulations of 
induced superselection rules are discussed in Sect. 5. 

2 Structure of the state space 

We start with a few mathematical notations. If 7i is a separable Hilbert space 
we use the following spaces of linear operators defined on Ti. 

B{l-L): The R-hnear space of all bounded selfadjoint operators A. The norm 
of this space is the operator norm 

TiTi): The R-linear space of all selfadjoint traceclass operators A. These 
operators have a pure point spectrum G R, i=l,2,..., with J2i \c(i\ < oo. The 
natural norm of this space is the trace norm ||v4||i = tryA+Z = X]j 
other norm, used in the following sections, is the Hilbert-Schmidt norm \\A\\2 — 
VtrA+A. These norms satisfy the inequalities ||v4|| < ||A||2 < jj^ljji. 

V{7i): The set of all statistical operators, i.e. positive traceclass operators 
W with a normalized trace, iiW = 1. 

V{7i): The set of all rank one projection operators P^. 
These sets satisfy the obvious inclusions V{n) C V{n) C T(7i) c B{n). 

Any state of a quantum system is represented by a statistical operator W e 
I){T-C), the elements of P(7i) thereby correspond to the pure states. Any observ- 
able is represented by an operator A G BiTi). The use of bounded operators 
is only an apparent restriction, since the whole information of an unbounded 
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observable can be recovered from its spectral resolution involving only bounded 
operators. The expectation of the observable A in the state W is the usual trace 
trWA, which is always finite for W e V{H) and A E B{H). Without additional 
knowledge about the structure of the system we have to assume that the set of 
all states corresponds exactly to T>(TC), and the set of all (bounded) observables 
to BiH). The state space V(TC) has an essential property: it is a convex set, i.e. 
Wi,W2e V{H) implies AiK^i + A2W2 e V{n) if Ai,2 > and Ai + A2 = 1. Any 
statistical operator W G ViTi) can be decomposed into pure states 

W = Y. ^nPl (1) 
n 

with G ViTi) and probabilities Wn > 0, J2n''^n = 1- An explicit example is 
the spectral decomposition of W. But there are many other possibilities, and we 
shall investigate that aspect in more detail. It is exactly this arbitrariness that 
does not allow a classical interpretation of quantum probability. 

Before we continue with the investigation of the quantum mechanical state 
space we introduce some notations concerning convex sets. A bounded and closed 
subset A4 of an R-linear space is convex if with two points xi, X2 E Ai also the 
connecting line segment Xxi + (1 — A)a;2, A G [0,1], belongs to M.. Those points of 

which do not lie on a line connecting two others, i.e. for which x = \xi + \x2 
with xi^2 G TVI is only possible if xi = X2 = x, are called extremal. The extremal 
points are always boundary points, but not necessarily all boundary points are 
extremal. The boundary of A4 will be denoted by dA4, the set of extremal points 
- the extremal boundary - by deM.- An important statement about convex sets is: 
Any point x E A4 can be represented by an integral over the extremal boundary 

x = [ ydX{y) (2) 

where dX{y) is a probability measure, i.e. a non-negative measure with / dX{y) = 
1, concentrated on (the closure of) dgAi. This representation has been derived 
in a rather general context by Choquet, see e.g. [11,12]. In the case of point 
measures the integral (^ yields the sum 

i 

with yi G deM. and probabilities < Aj < = 1- The representations 

(0) and(^ are only unique if is a simplex. An n-dimensional simplex Ai 
is the closed convex set generated by n+1 points yi,i = l,...,n + 1, where 
the connecting lines yi — yn+i,i = l,.-.,n are linearly independent. At ={x = 
J2i=i,...,n^iyi\K ^ 0,X]iAj = 1}. In two dimensions this is a triangle, in three 
dimensions a tetrahedron. The extremal boundary of a simplex is the set of its 
vertices, and the weights Aj are the uniquely defined barycentric coordinates of 
the point x E Ai with respect to the vertices yi. Another class of convex sets 
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are balls = {x G R", |x| < 1}. It is easily seen that in this case the extreme 
boundary coincides with the boundary S*""^ = {x G R", |x| = 1}. For balls the 
representations (Q) or are highly non-unique, since any y G deM. can show up 
on the right hand side with positive weight. 

The essential difference between the state space of classical mechanics and 
that of quantum mechanics is: 

- The state space of classical mechanics is an (infinite dimensional) simplex with 
a unique representation (|^). 

- The state space of quantum mechanics is not a simplex, its structure is closer 
to that of a ball. The decompositions (H) or (|) into extremal elements, i.e. pure 
states as explained below, are highly non-unique. 

The state space of classical mechanics is the set of all probability measures 
on the phase space S C R" of a given system. The extremal boundary 
of that set is given by all Dirac measures {5r?(0^"^ = ~ vi)d"'^\r] G S}. 
The representation (0) is the identity d^{^) = J 6rj{^)d^{ri), which leads again 
to the originally given measure. For the mathematical proof that the space of 
probability measures is an infinite dimensional simplex with a unique Choquet 
representation see [11,12]. 

We have already seen that the state space of quantum mechanics is convex. 
The boundary of this set is formed by those statistical operators which have 
at least one eigenvalue zero. The condition for extremal elements, i.e. W = 
\Wi = 1 14^2 has to imply Wi = W2, is met by all rank one projection operators 
W G V(H). The decomposition (|T]) is therefore the Choquet representation of 
an element of a convex set. As a consequence of the superposition principle this 
representation is highly non-unique, see e.g. Sect. 2. 3 of [13]. As explicit example 
we consider the state space of a spin-| system, i.e. V{C'^). Since all statistical 
operators on with one eigenvalue zero are rank one projectors, the boundary 
of the state space coincides with the extreme boundary. More explicitly, any 
operator p G P(C^) can be represented with a real polarization vector 'jf in the 
unit ball = {jf G R^||"p^| < 1} as 

Pit) ^ (1 + ^t) (4) 

where 1 is the unit 2x2 matrix and "o^ = (cri, a2, C3) are the Pauli matices. More- 
over, a simple calculation leads to the identity ||p(l?i) — ^("^^2)111 = ["P^^i ~ 1^2l- 
Hence the topology of P(C^) given by the trace norm agrees with the euclidean 
metric of the polarization vectors. The set P(C^) is therefore isomorphic to 
the unit ball B^. The extremal elements corresponding to 7^ G S"^, are exactly 
the projection operators of the representation (|I]). Since Xipi'P'i) + ^2P{1^2) = 
p{^il^i + ^21^2) if -^1,2 > 0, A1 + A2 = 1, the representation (|l]) corresponds there- 
fore to the Choquet representation (|^) of the unit ball B^. The non-uniqueness 
of that representation is therefore obvious. 
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In higher dimensions the geometrical picture of V{T-C), dirsiH > 3, is more 
comphcated. The boundary is much larger than the extremal boundary, but the 
arbitrariness of the representation (|l]) remains. To see this one can choose any 
two dimensional subspace of the spectral decomposition of W with at least one 
non-zero eigenvalue. Then this part of the spectral representation corresponds 
up to normalization to P(C^), and we can modify it with all the arbitrariness 
seen above. 

3 Superselection rules 

The arbitrariness of the decomposition (Q) originates in the superposition princi- 
ple. In quantum mechanics, especially in quantum field theory, the superposition 
principle can be restricted by superselection rules. Here we cannot discuss the 
arguments to establish such rules, for that purpose see e.g. [5,6], and also Chap. 6 
of [10], or to refute them, see e.g. [14]. Here we only investigate the consequences 
for the structure of the state space. In a theory with superselection rules like the 
charge superselection rule, the Hilbert space H splits into orthogonal superselec- 
tion sectors Hm, m G M, such that H = Q)m 'Hm- Pure states with charge m (in 
appropriate normalization) are then represented by vectors in Tim-, and super- 
positions of vectors with different charges have no physical interpretation. The 
projection operators onto the orthogonal subspaces Tim satisfy PmPn = ^mn 
and X)m Pm = I- The set of states is reduced to those statistical operators which 
satisfy PmW = WP^ for all projection operators Pm, rn G M. The state space of 
the system is then = {W e V{n)\WPm = P^W, m e M}, and all statistical 
operators satisfy the identity 

W = J2PmWPm. (5) 
m 

A mathematical equivalent statement is that all observables of such a theory 
commute with the projection operators Pm,m G M. Superselection rules of this 
type will be called "kinematical superselection rules" to contrast them against 
the "dynamically induced superselection rules" to be discussed in the following 
section. 

As a consequence of (|^) the extremal boundary of T)^ decomposes into = 
UmV^ with = G P{n)\P^n C Hm}- The decomposition of a statis- 
tical operator W G V''^ into pure states now reads (here the sum over i can be 
substituted by an integral over the set P^) W = J2m,i ^m,iPm,i ^^^^ Pm,i ^ 
and probabilities A^.i > 0, J2m,i ^m,i = 1- Since the representation of an element 
of ViTirn) by pure states is not uniquely given (if dimHm > 2), the right hand 
side of of this decomposition is again highly non-unique. But as a consequence 
of the structure of the state space the total probability with respect to any 
of the subsets P^, i.e. J2i ^m,i, has no ambiguity since J2i ^m,i = trl^P^, where 
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the right hand side is uniquely defined. Hence superselection rules allow to speak 
about objective "properties" of a quantum system. These properties show up 
with the probabilities trWPm irrespective of what other specifications (measure- 
ments) are made. In the language of observables these properties are represented 
by the self adjoint operators Pm, which commute with all observables of the sys- 
tem. The importance of superselection rules for the transition from quantum 
probability to classical probability is obvious. But there remains an essential 
problem: Only very few superselection rules can be found in quantum mechanics 
that are compatible with the mathematical structure and with experiment. A 
satisfactory solution to this problem is the emergence of effective superselection 
rules induced by the interaction with the environment. 

4 Dynamics of subsystems 

In the following we consider an "open system" , i.e. a system 5" which interacts 
with an "environment" E, such that the total system S + E satisfies the usual 
Hamiltonian dynamics. The system S is singled out by the fact that all observa- 
tions refer only to this subsystem. The Hilbert space TCs+e of the total system 

5 + E is the tensor space TCs ® 7i_E of the Hilbert spaces for S and for E. We 
assume that the only observables at our disposal are the operators Ie where 
A G BiJ-is) is an arbitrary bounded selfadjoint operator on Tis- If the state of the 
total system is G V^TCs+e), then all expectation values trVF(y4 ® Ie) can be 
calculated from the reduced statistical operator p = ti eW which is an element of 
V{T-Cs) defined such that the expectation values satisfy tr^p^ = tTs+EW{A(^lE)- 
Since all information about a physical subsystem is given by a statistical opera- 
tor, we shall here refer to the statistical operator as the "state" of the subsystem. 
(The state of the total system cannot be recovered from these "states" of its 
subsystems.) 

As mentioned above we assume the usual Hamiltonian dynamics for the total 
system, i.e. W{t) = U{t)W{0)U~^{t) with the unitary group U{t), generated by 
the total Hamiltonian. Except for the trivial case that S and E do not interact, 
the dynamics of the reduced statistical operator p{t) = tTEU{t)W{0)U~^{t) is no 
longer unitary. It is the purpose of the second part of this talk to evaluate this 
dynamics in some detail. The essential result is that this dynamics can produce 
effective superselection sectors, i.e. 

p{t)=J2PmPit)Pm (6) 

m 

in sufficiently short time with a set of projection operators Pm-,^ ^ M, which 
correspond to a superselection structure (|^). In Sect.|| we shall give more precise 
formulations of (^. The suppression of the off-diagonal terms of the statistical 
operator in (^ is essential for the emergence of classical properties in quantum 
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mechanics. This suppression can never be understood by semiclassical approxi- 
mations alone. 



5 Solvable models 



The statement is so far rather vague since it does not specify the asymptotics. 
The following examples show what type of asymptotics is possible in principle. 
For a restricted class of models an estimate 



\\p{t)-Y.p^p{t)p^ 



< c^{i + 5\t\y 



(7) 



with the Hilbert-Schmidt norm is possible for arbitrary p(0) = p E V{T-ls). 
The constants 7 > 0, (5 > and > do not depend on p. Since one can 
achieve large values of 7 and/or small values of the constant C-y, these dynamically 
induced superselection sectors PnTi. cannot be distinguished practically from the 
kinematical superselection sectors @. For a subclass of these models one can even 
derive an estimate of the type (|^) with the stronger trace norm. A consequence of 
(1^ is that the transition between a superselection sector P„7^ and its complement 



PnTi. := (/ — Pn)^^. is suppressed by 



< C^(l + \t\) "I again with the 



Pnp{t)Pn 

Hilbert-Schmidt norm. If P is a projection operator with finite rank on a 
subspace of i.e. F„F = FP„ = F, then Pnp{t)P is a traceclass operator 



with Pnp{t)P ^ < ^/NC^{l + \t\)~''' , and the transition between PTi and all other 
sectors Pm'H, m 7^ n, is uniformly suppressed in trace norm. The investigation 
of the models also shows that unfortunately such simple estimates are rather 
unstable against slight modifications of the models. 

The models have the following structure. The Hilbert space is Hs+e = 'Hs ® 
T-Ce- The total Hamiltonian has the form 



H = Hs® Ie + Is® He + Vs(S>Ve + V 



(8) 



where Hs is the Hamiltonian of S, He is the Hamiltonian of E, V5 (S> Ve is the 
interaction term between S and E with selfadjoint operators on Tis and Ve 
on TIei aiid y is a possible additional scattering potential. For all models we 
assume that 



1) The operators He and Ve commute, [He, Ve] = O, hence [Is 
O. 

2) The operator V5 has a pure point spectrum 

^5 = 51 ^rnPm 



He, Vs 



Vk 



(9) 



where the projection operators form a complete set, PmPn = ^m,nPn, and 
Hm P'tn = P We assume that the eigenvalues are separated with the lower 
bound I Am — An I > 5 > if m 7^ n. 
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3) The operator Ve has an (absolutely) continuous spectrum. 

Remark. The continuous spectrum of Ve is needed to obtain simple estimates 
for if: ^ oo. But one could also allow an operator with point spectrum (as done 
in [8]), if the spacing of the eigenvalues is sufficiently small. Then the norm in 
(1^ is an almost periodic function and the inequality is only correct for a finite 
time interval <t <T. But T can be large enough for all practical purposes. 



5.1 The Araki-Zurek model 

The first solvable models to discuss the reduced dynamics have been given by 
Araki [7] and Zurek [8], and the following construction is essentially based on 
these papers. In addition to the specifications made above, we demand: 

4) The operators Hs and Vs commute, the potential V vanishes, i.e. [Hs, Vs] = 
= V. 

For an originally factorizing state W = p®uj we calculate with U{t) = exp{~iHt) 

Pit) = tTEU{t)WU'-{t) = e-^^^* Yl PmPPne'^^'^XmAt) (10) 

m,n 

with Xm,ni^) = (e-*(^™-^")^^^*u;) , see e.g. Sect. 7.6 of [10]. The trace XmA^) 
vanishes for |t| ^ cxd if m 7^ n since Ve has an absolutely continuous spectrum. 
Under additional smoothness assumptions on u (which do not restrict p G V(Ti.s)) 
we derive for x(^) := tr (^e~*^^*a;^ the estimate \x{i)\ < + \t\)~'^. Here 7 can 
be arbitrarily large if a; is a sufficiently differentiable function in the spectral 
representation of Ve (and vanishes at the boundary points of the spectrum). 
This estimate leads to the upper bound 

\XmAt)\<C,{l+S\t\)-^ (11) 

if I Am — A„| > S > 0, and we obtain the estimate (0). If we have only a finite 
number of eigenvalues in (|^), a bound of the type (|^) with the stronger trace 
norm can be derived. 

This result depends on the reference state u only via the decrease of x(t). 
We could have chosen a more general initial state W = J^^ Pfi® ^ 'P'iJ^s+E) 
with p^ G T(Hg), G V(He)- The operators p^ need not to be positive 
separately, but we have of course p = Y^^Pti ^ T^O^s)- Then (^ is still valid if 
Yl,^ tr |^e~*'-'^'"~^"''^^*to'^^ satisfies a uniform estimate (pHI) . Hence the emergence 
of the superselection sectors Pn'Hs is not sensitive to the initial conditions. 



5.2 A spin model 

So far we have made the rather restrictive assumption V = [Hs,Vs] = O such 
that the interaction commutes with the free Hamiltonian Hs ®Ie- The following 
model has still V = but [i/5, V^] does no longer vanish. 
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The Hilbert space of the model is Hs+e = 'Hs ® 'He with Tis = and 
TIe = C?{^)- The Hamiltonian has the form (j^) with V = O and the following 
specifications 

Hsi^ = C(i'^)i^, with "(f G R^, Pauh matrices, i/j E C"^, 
HEf{x) = hx^ f{x) with a positive constant 6 > 0, f{x) G £^(R), 
Vs = Xcr-i, with a real coupling parameter A, 
VEfix) = xf{x),fix)eC'{R). 

The statistical operator of the spin-| system is a spin density matrix (^) with 
a polarization vector 'p' E B^. For the total system we assume an initial state 
W = piJ^) where a; is a statistical operator on >C^(R) with a smooth integral 
kernel uj{x,y). The time evolution U{t) = exp{—iHt) with the total Hamiltonian 
(H) then leads to the state U{t)WU^{t). The diagonal part of the integral kernel of 
this statistical operator is {U{t)WU'^{t)) {x,x) = u^x, x)p {Rx{t)'jf) . Here Rx{t) 
is the rotation induced by the SU{2) matrix part of U{t), i.e. exp{—iht), with 
h = ('a*' a^) + Xxa^. The reduced statistical operator p(t) of the spin-i system 
is calculated as pit) = tTEU{t)WU~^{t) = J dxu{x,x)p {Rx(t)jf'). For the initial 
state p = tr^iy = p{'jf) the reduced dynamics then leads to a statistical operator 
p(t), which asymptotically approaches p{lf), where the polarization vector if is 
given by the linear mapping 

if = MJ)' := J dxuj{x,x)li{x) {jflt{x)) (12) 

with the axis lt{x) of the rotation Rx{t). Under appropriate conditions for the 
initial state uj of the environment, the difference p{t) — p{lf) can be uniformly 
estimated by \\pit) —pi~f')\\i < c(l + Itl)""^. The mapping (|12]) is a symmetric 
contraction on R'^. We can distinguish two cases: 

1) If II 63 the mapping ([12|) reduces to Mjf = 63 {e^jf), and we obtain the 
results discussed in Sect.^.l|. The condition |M"jf | = is satisfied for "jf || 63, 
and only in this case p{t) is not affected by the decoherence. 

2) If 1^ has components orthogonal to 63, also the direction of depends on 

and lAfp^l < 1 1?^| holds for all vectors ^ . 

In the second case there are no projection operators which commute with all 
operators p^Mj^), Jf G B^. A superselection rule of the type (^) can hold in 
some approximative sense only if a\ + a\ <^ a\. In that case has very 

small components orthogonal to 63, and the off-diagonal matrix elements of the 
operators p{MJ^), G B^, are negligible. 

5.3 Models with scattering 

We assume again a Hamiltonian (^) with [Hs, Vs] = O as in the Araki-Zurek 
model, but now with an additional scattering potential V defined on the full 
Hilbert space Hs+e = 'Hs ®'He- Under appropriate conditions on V the wave 
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operator exists as strong limit Q = lim(_>oo e*''^*e~*^"* on Tis+E-, where Hq is the 
Hamiltonian of the Araki-Zurek model, = Hs ® Ie + Is ® He + Vs® Ve- (We 
assume for simplicity that there are no bound states such that Q'^ = Q^^.) Then 
the time evolution U{t) = exp{—iHt) behaves asymptotically as Uo{t)Q'^ with 
[/o(t) = exp^—iHot) . More precisely we have for all W G 



lim 

t— >oo 



u{t)wu+{t) -Uo{t)n+wnu^{t) =o (i3) 



in the trace norm. For initial states W = p(0) ® uj with smooth u, and for 
sufficiently regular scattering potentials the statistical operator Q~^WQ has only 
smooth contributions in T(H), and the reduced trace ti EUo{t)Q'^WQUQ (t) yields 
the induced superselection sectors Pm'Hs of Sect.pTl|. But then (|13|) implies for 
p{t) = tTEU{t)WU+{t) the asymptotics hmt^^\\p{t) - En, PmPit)Pm\\2 = 0. 
Hence p{t) has again the induced superselection sectors which originate from 
the spectrum of Vs. But in contrast to (0) this bound is not uniform in the 
initial state p(0), since scattering does not allow a uniform estimate for (plSf) . 



5.4 Concluding remarks 

The investigation of the models proves that the uniform emergence (0) of supers- 
election sectors is consistent with the mathematical rules of quantum mechanics. 
But this result depends on rather restrictive assumptions on the Hamiltonian. For 
more realistic models we can nevertheless expect a strong quantitative suppres- 
sion of the off-diagonal elements of the statistical operator, as already proposed 
in [8]. 
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